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O ABSTRACT 

c . 

d . The stratification-driven Hall instability in a weakly ionized polytropic plasma is 

investigated in the local approximation within an equilibrium Keplerian disk of a small 
aspect ratio e. The leading order of the asymptotic expansions in e is applied to both 
equilibrium and perturbation problems. The equilibrium disk with an embedded purely 
qq ■ toroidal magnetic field is found to be stable to radial, and unstable to vertical short-wave 

perturbations. The marginal stability surface was found in the space of the local Hall 
and inverse plasma beta parameters, as well as the free parameter of the model related 
to the total current. To estimate the minimal values of the equilibrium magnetic field 
that leads to the instability, it was constructed as a sum of a current free magnetic field 
and the simplest approximation for magnetic field created by the distributed electric 
Qh' current. 
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1. Introduction 

The dynamics of thin rotating gaseous disks under the influence of a magnetic field is of great 
importance to numerous astrophysical phenomena. Of particular importance is the study of the 
various instabilities that may arise within the disks whose inverse growth rates are comparable 
to the rotation period. The rediscovery of the magneto-rotational instability (MRI) within astro- 
physical context (Balbus and Hawley 1991) gave the sign to a series of extensive investigations of 
hydromagnetic instabilities of rotating disks. 

While most of the research has focused on the magnetohydro dynamics (MHD) description of 
MRIs, the importance of the Hall electric field to such astrophysical objects as protoplanetary disks 
has been pointed out by Wardle (1999). Since then, such works as Balbus and Terquem (2001), 
Salmeron and Wardle (2003, 2005), Desch (2004), Urpin and Rudiger (2005), and Rudiger and 
Kitchatinov (2005) have shed light on the role of the Hall effect in the modification of the MRIs. 
However, in addition to modifying the MRIs, the Hall electromotive force gives rise to a new family 
of instabilities in a density-stratified environment. That new family of instabilities is characterized 
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by the Hall parameter £ of the order of unit (£ is the ratio of the Hall drift velocity, Vhd, to 
the characteristic velocity, where Vhd = If^i/H*, li and fij are the inertial length and Larmor 
frequency of the ions, and H* is the density inhomogeneity length, see Huba (1991), Liverts and 
Mond (2004), Shtemler and Mond (2006)). As shown by Liverts and Mond (2004) and by Kolberg 
et al. (2005), the Hall electric field combined with density stratification, gives rise to two modes of 
wave propagation. The first one is the stable fast magnetic penetration mode, while the second is 
a slow quasi-electrostatic mode that may become unstable for short enough density inhomogeneity 
length H*. The latter is termed the Hall instability, and the aim of the present work is to introduce 
it within astrophysical context. 

Although the magnetic field configuration in real disks is largely unknown, observations and 
numerical simulations indicate that the toroidal magnetic component may be of the same order 
of magnitude as, or some times dominate the poloidal field due to the differential rotation of the 
disk, even if the magnetic field is created continuously from an initial purely poloidal field (see 
e.g. Terquem and Papaloizou 1996, Papaloizou and Terquem 1997, Hawley and Krolik 2002, Proga 
2003). 

The present study of the Hall instability is carried out for thin equilibrium Keplerian disks 
with embedded purely toroidal magnetic fields. While the accepted approach to the steady state 
of the rotating disks in most astrophysics-related hydromagnetic stability researches is that of a 
"cylindrical disk", a more realistic model of thin disks is adopted in the current work. Thus, 
an asymptotic approach developed for equilibrium non-magnetized rotating disks (Regev, 1983, 
Kluzniak and Kita 2000, Umurhan et al 2006), is employed here in order to construct a family 
of steady-state solutions for the rotating disks in a toroidal magnetic field within the Hall-MHD 
model, and to their further study of their stability. It is shown that in such radially stratified 
systems the Hall instability indeed occurs with the inverse growth rate of the order of the rotation 
period. 

The paper is organized as follows. The dimensionless governing equations are presented in the 
next Section. In Section 3 the Hall-MHD equilibrium configuration of a rotating disk subjected to a 
gravitational potential of a central body and a toroidal magnetic field, is described in the thin-disk 
limit. Linear analysis of the Hall stability of toroidal magnetic fields in thin equilibrium Keplerian 
disks and the results of calculations for both arbitrary and specific magnetic configurations are 
presented in Section 4. Summary and conclusions are given in Section 5. 

2. The physical model. Basic equations 

2.1. Dimensional basic equations 

The stability of radially and axially stratified thin rotating disks threaded by a toroidal mag- 
netic field under short-wave perturbations is considered. Viscosity and radiation effects are ignored. 
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The Hall electric field is taken into account in the generalized Ohm's law which is derived from the 
momentum equation for the electrons fluid neglecting the electron inertia and pressure. Under the 
assumptions mentioned above, the dynamical equations that describe the partially ionized plasmas 
are: 

_ + V-(nV) = 0, (la) 
_DV 1 

rriin—— = —VP + - j x B — m,jnV$, (lb) 
Dt c 

§- t (Pn^) = 0, (lc) 

<9B 

— + cV x E = 0, V • B = 0, (Id) 

E = --V x B + llii^, j = x B. (le) 
c ec n e An 

Cylindrical coordinates {r, 6, z} are adopted throughout the paper with the associated unit basic 
vectors {i r , ig, i z }; t is time; <& is the gravitational potential of the central object; &(R) = —GM/R; 
R 2 = r 2 + z 2 ; G is the gravitational constant; and M is the total mass of the central object; c is 
the speed of light; e is the electron charge; E and B are the Hall electric and magnetic fields; j 
is the electrical current density; V is the plasma velocity; D/Dt = d/dt + (V • V) is the material 
derivative; n is the number density; P = P e + Pj + P n is the total plasma pressure; P^ and are 
the species pressures and masses (k = e,i,n); subscripts e, i and n denote the electrons, ions and 
neutrals. In compliance with the observation data for protoplanetary disks, the density decreases 
outward the center of the d isk according to the law n ~ \ jr a with the density exponent a ~ 2.75 



(jGuilloteau Dutrev 19981 ) . In a standard polytropic model of plasmas a = 3/(7 — 1), this yields 
the specific heat ratio 7 ~ 2.1 that is somewhat higher than the standard value 7 = 5/3 used here. 
Since the plasma is assumed to be quasi-neutral and weakly ionized 

n e & rii i=a an n , n n n . (2) 

This indicates the increasing role of the Hall effect in the generalized Ohm law with vanishing 
ionization degree, a = n e j(n e + n n ) « 1. 



2.2. Scaling procedure 

The physical variables are now transformed into a dimensionless form: 

fnd = f/f*, (3) 

where / and f n d stand for any physical dimensional and non-dimensional variables, while the 
characteristic scales /* are defined as follows: 

V* = r2*r*, t* = — , = V 2 , m* = mi, n* = n n , (4a) 

1 ijj. 
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P* = K(m*n*) 7 , j* = — - , E* = Yl^_ (4b) 

47T c 

Here $7* = (GM/r 3 ) 1 ' 2 is the value of the Keplerian angular velocity of fluid at the characteristic 
radius r*; K is the dimensional constant in the steady-state dimensional polytropic law. The 
characteristic values of the magnetic field B* and radius r* are the dimensional parameters of 
the problem. This yields the following dimensionless system (omitting the subscript nd in non- 
dimensional variables) : 

^ + V • (nV) = 0, (5a) 
DV 11 1 

"In =-Ml yp + Wl ixB - nV "- ^^ = -J^T7W (5b) 

= »• (5c) 

— + VxE = 0, V • B = 0, (5d) 
at 

E = -VxB + ?— , j = VxB, (5e) 
n 

where Ms and (3 are the Mach number and plasma beta, and £ is the Hall coefficient: 

M S = — P = 4tt— £= ni ( li Y= B *° (6) 



cs* = \J 7-P* I n* , h = c/ujpi and Qi = eB*/(rriic) are the inertial length and the Larmor frequency 
of ions, respectively, u) p i = y / 4Tre 2 rii/mi is the plasma frequency of the ions. 

A common property of thin Keplerian disks is their highly compressible motion with large 
Mach numbers Ms 

_L = C «1, e = El : (7) 
M s n 

where e is the aspect ratio of the disk; is a characteristic dimensional thickness of the disk. It 
is convenient therefore to use a stretched vertical coordinate 

z H 

C = -,h = tL. (8) 



3. Equilibrium solution in thin-disk limit for toroidal magnetic configuration 

The steady-state equilibrium of rotating magnetized disks may now be obtained to leading 
order in e by writing all physical quantities as asymptotic expansions in small e (similar to Regev, 
1983, Kluzniak and Kita 2000, Umurhan et al 2006 in their analysis of accretion disks, see also 
references therein). This yields for the gravitational potential in the Keplerian portion of the disk: 

Hr,0 = - 1 - + e^ 3+ O(e% (9) 
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To zeroth order in e , the toroidal velocity Vg is described by the Keplerian law: 

V e = fi(r)r, n{r) = r~ 3/2 , (10) 

the rest of the velocity components are of higher order in e, and their input into the stability 
analysis may thus be neglected. Let us assume now that the magnetic field is purely toroidal and 
depends on both r and z: 

B = B e (r,z)i e , j = VxB, (11a) 

. d 1. d . ,. n . , r laSfl IdirBg). 
V = i r — + -i 2 — , j = {j r ,0,j z } = {— " 0,-^— ^}. lib 

Or e o( e o( r Or 

Since the vector product of the two toroidal vectors V and B is zero, it follows from system 
(1) that the Hall electric field has the potential 0: 

j x B 



n 



= V0 (12a) 



or, equivalently, 

d4> Bg djrBg) d<t> dBg 
n Tr = ~ V^-' n d( = ~ ~d(~' (12b) 

Substituting Eqs. (12) into the vertical component of the momentum equation (5b) yields in 
the leading order in e: 

i i 



P 

cs = A/7 — j P = ri y , n = u 
n 



r,() , h 2 -( 2 



+ 



7-1 / ry _ \ \ 7-1 



'-IV 1 (13> 



(3 ' r 3 

Here h is the disk thickness, where n and <p equal zero 

n = 0, = at C = h, (14) 
h(r) is scaled so that h(l) = 1; /i(r) should be specified to close the problem (see Section 4). 

Finally, a general expression for the toroidal magnetic field as well as for the Hall potential 
may be derived by solving the induction equation that assumes the following form: 

d{r 2 n) d{r 2 B 2 e ) d(r 2 n) d{r 2 B 2 e ) 

~8X d~r o~r lK~ ' ( ' 

This equation is satisfied for rBg that is an arbitrary function of N (Kadomtzev 1976), where N 
has the meaning of the inertial moment density (see additionally Section 4.3): 

B e = r~ 1 I{N), N = r 2 n. (16) 

This yields for the Hall potential 

m) _£mm dN , (17) 
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where the arbitrary constant in the potential is chosen from the condition that the potential is zero 
at the disk edge where n = 0, i.e. (f)(0) = 0; the upper dot denotes a derivative with respect to 
the argument. If the arbitrary function I(N) is given, then Eqs. (13) and (17) constitute implicit 
algebraic equations for the number density n and the potential 4>. 



4. Hall instability of thin Keplerian disks 

The stability properties of the axially symmetric steady-state equilibrium solutions for thin 
Keplerian disks are investigated now under small axially symmetric perturbations, and find the 
dispersion relation in the short wave limit. 



4.1. Thin disk approximation for perturbed problem 

We start by linearizing the Hall MHD equations (5) about the equilibrium solution within the 
Keplerian portion of the thin disk. For that purpose, the perturbed variables are given by: 

F(r,z,t)=F(r,()+F'(r,(,t). (18) 

Here F stands for any of the physical variables, F denotes the equilibrium value, while F' denotes 
the perturbations of the following form: 

B'(r,z)=B' e i e , j'« {1*3,0, 23}, V' = {Q,0X}. (19) 
Inserting the relations (18)-(19) into Eqs. (5), linearizing them about the equilibrium s olution 



and applying the principle of the least possible degeneracy of the problem (jVan Dyke 1964 I ) to the 
result, we obtain that B' e is of the order e°, while the axial perturbed velocity V z is of the order 
e 1 and should be rescaled. Furthermore, we rescale the Hall parameter, keeping in mind that it 
should be of the order of e , since its higher and lower limits in e lead to more degenerate problems. 
Thus, rescaling the axial perturbed velocity and Hall parameter, we derive 

y „ f Q./2 

v - = T' « = ! s nTTT (20> 

Here the rescaled Hall parameter £ is equal to the ratio of the Hall drift velocity, Vhd, to the 
characteristic velocity, V* = Q*r*, and the characteristic disk thickness H* = er* (see Eq. (7)) is 
the density inhomogeneity length. It is noted that the limit £ << 1 yields the MHD model. Finally 
this leads to the following linearized system in the leading order in e: 



(21a) 
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8V: _ dP> ldB e 1 dB' e 

n ^-~^Q- (3^ Be ~ (21b) 
dP' 2 dn' 

^- C % = °< < 21c > 

dB' e d{B' d V>) ;B e \l d(rB e ) dn> 1 d(r 2 n) dB' e , dndB' e dB e dn>\ ~ 



dt dQ n 2 \ r dr d( r 2 dr d( d( dr dQ dr 

Here the perturbed toroidal magnetic field identically satisfies the relation V • B' = 0; ^ denotes 
the terms in the magnetic induction equation that do not include derivatives of perturbations. The 
latter will be neglected below in the short wave limit, and it is not presented here for brevity. 



4.2. Dispersion relation in the short wave limit 

The perturbed quantities are assumed to satisfy the following conditions: 
, 1 dF' , , 1 dF' , , 1 dF' , 1 dF 1 dF 1 

| |^| | rsj — ^ 1 . (^^/ 



F' dt F' d( F' dr F d( F dr r 

This allows us to assume the following local approximation for the perturbations: 

F'(r,C,t) = f exp(-ioot + ik r r + ik c (). (23) 

Here /' denotes the amplitude of fluctuations about the equilibrium; oo is the (generally complex) 
frequency of perturbations; k = {k r ,k^} is the wave vector in the scaled coordinates {r, £}. It 
should be distinguished from the wave vector k = {k r ,k z } in the {r, z} coordinates, k z = k^/e. 
Condition (22) may now be written as 

, _ . 1 OF . . 1 dF . 1 1 
w ~ ~ kr > I H j-^ |~ - ~ 1. (24) 

Using relations (23)- (24), the system (21) may be presented as the following system of homo- 
geneous linear algebraic equations: 

- Cri + nV' z = 0, (25a) 
- CnV' z = -P' - p- x B e B' e , v; = 0, (25b) 
P' = c 2 s n', (25c) 

Here the scaled phase velocity C = uo/k^ in the scaled coordinates {r, (} is related to the phase 
velocity C in the cylindrical coordinates {r, z} as C = C/e, and 

<4 = <4r - <4 l - = T k= T k= <2a> 
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Note that since the radial wave number k r is contained only in the Hall term in the magnetic 
induction Eq. (25d), in the leading order in e the standard MHD problem is free from k r . 

The system of equations (25) has non-trivial solutions if the following cubic eigenvalue equation 
for the local value of the rescaled phase velocity c is satisfied: 

Xc 3 +Xc 2 -c-q = 0. (27) 

Parameters X and q are expressed through the local values of the plasma beta b, Hall coefficient x 
and new parameter g whose the physical meaning will be elucidated later 

c = , X = x , q= —r— x , 28a) 

xcs 1 + 1 1 + 1 

6-.^, ,=iA4f. («•) 

c s n csn z r z a/3 a LnJy 

where Bq and N are determined in Eqs. (16). Straightforward analysis of dispersion relation (27) 
(such dispersion relation for the Hall instabili ty of plasma in slab geometry has been derived and 



investigated in (jBrushlinskii &: Morozov 19801 )) reveals that the system is stable for < q < 1, 
while the instability may occur for the following two regimes: 

1) q > 1 and X~ < X < X + or 2) q < and X > X~ . (29) 

Here the first instability regime corresponds to g > 1, and the second one to g < —b (see Eqs. 
(28a)); X~ and X + are the roots of the following quadratic equation: 

4qX 2 + (1 + 18q - 27q 2 )X + 4 = 0. (30) 

Although Eq. (30) provides for the dependence of the local Hall parameter X vs q, it is more 
convenient to use X as a function of the local plasma beta b and the free parameter g, which may 
be related to the natural parameter of the disk - total current through the disk (see Subsection 
4.3). As a first result, it is noticed that the current free magnetic configuration with g = is stable 
since < q < 1. In addition, it is noted that the quantities X and q (or, equivalently, X, g and b) 
which determine the stability properties of the disk, depend on the local values of the equilibrium 
number density and magnetic field which are determined in Section 3 up to the free function I(N) 
of the inertial moment density. 

Figure 1 depicts the family of the marginal stability curves which correspond to the real roots 
of Eq. (30). These curves are presented in the plane X — b^ 1 of the local Hall parameter and inverse 
plasma beta with a single parameter of the family, g. According to the relations (29), there are two 
instability regimes (Figs, la and lb). Thus, in Fig. la the values of X located between the upper 
and bottom branches of the marginal stability curves give rise to the first regime of instability, 
while the disk is stable if X falls outside that interval. In the limit 6 _1 — > 0, we have q — > 1, and 
Eq. (30) has a double root X = 1. This means that in the limit of high plasma beta the instability 
interval becomes very narrow and the disk is stable. In Fig. lb, the values of X located higher the 
line X~ correspond to the second instability regime in (29) (q < mode), and the disk is stable 
otherwise. 
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4.3. Model equilibrium solutions for toroidal magnetic fields 



The marginal stability curves in Fig. 1 provide for only a qualitative description of the disk 
stability, since they are depicted in terms of the local parameters which are determined up to the 
free function I(N). To estimate qualitatively the stability characteristics of the system a specific 
example of a toroidal magnetic field is considered. The free function I(N) is presented as a sum 
of the current free magnetic field (created by a current localized along the disk axis), and the 
simplest quadratic approximation for the magnetic field created by an electric current distributed 
over Keplerian portion of the disk: 



I(N) =I + I: 



N 2 



2 > 



N 



r n, 



(31) 



where Io, I2 are constant coefficients; uq = n(l,0) is the value of number density at the reference 
point in the midplane {r = 1, z = 0} (analogous notations will be used below for the Hall electric 
potential 4>o = </>(l,0), free parameter go = g(l,Q), Hall coefficient £0 = £(1,0), etc). The absence 
of a linear term in N in Eq. (31) is a result of the requirement that the Hall electric potential <p is 
finite at the disk edge z = h, which in turn allows the number density to be zero. 

Substituting Eq. (31) into Eq. (17) and using the normalization condition Bg = 1 at the 
reference point {r = 1, z = 0} yield the following expression for the Hall electric potential: 



<f>(N) 



-21 



N 
2_ 2 



Io + 



3 n\ 



h + h 



1. 



(32) 



Definitions of g and I in Eqs. (28) and (31) relate the values of I2, go and the total current through 
the disk I t : 



50 = 2/ 2 , It = 2tt 



1 



1 



N 2 



11? 



go 
2 



(33) 



where = {r 2 n m ) r ^ c>:i is constant (see Figs. 2 and 3). This leaves go or, equivalently, the 
natural parameter of the disk - total current It as the single free parameter of the model. Equations 
(13) and (32)-(33) provide for the magnetic configuration in the equilibrium Keplerian disk. To 
overcome the presence of the arbitrary semi-thickness h{r) in relations (13), the small slope of the 
disk edge is assum ed (that is supported by the observation data for some protoplanetary disks 
(jCalvet et al. 20021 )). such that h(r) = 1. This yields at the midplane: 



n m (r) = v[f3 Vm + -?J , <t>m{r) 



7-1 



9or n m 



2 J 6 I no 



(34) 



Equations (34) form a coupled set of nonlinear equations for n m and <f> m , while no is determined 



from the following two nonlinear equations for no 

n = v \ P~ l <j>o + 1 



n. 



t (l) and 4> = <p m (l): 



1 

7-1 



go_ 

n 



go 
3 



(35) 
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Numerical solution of Eqs. (34) and (35) reveals two main types of equilibria that correspond 
to two possible signs of go, i.e. signs of the contributions of the distributed current to the total 
current. The midplane equilibrium magnetic field Bg m , number density n m , distributed electric 
current j zm and inertial moment density N m , are presented in Figs. 2 and 3 for positive and 
negative go, respectively. In both cases the equilibria are characterized by a distributed electric 
current localized at some radius that is moving outwards the disk center with increasing plasma 
beta. In addition, while for positive go both the equilibrium number density and magnetic field 
are decreasing monotonically to zero at a finite plasma beta, equilibria with negative go exhibit a 
maximum in the number density that corresponds to rings of denser material, as well as a maximum 
in the magnetic field profile. Finally, for go < the equilibrium solutions exist only for inverse 
plasma beta larger than a critical value corresponding to the vertical slope (/? _1 > 0.06 in Fig. 3), 
otherwise the equilibrium solution cannot be described by a single-valued function. Furthermore, 
Figs. 2d and 3d demonstrate that the inertial moment density N m tends to a non-zero constant at 
infinity through a local minimum (critical point) where d r N m =0. 



4.4. Hall instability of the model equilibrium 

Using the values of the free function I(N) in the magnetic field determined in the previous 
subsection, the family of the marginal stability curves presented in Fig. 1 may be recalculated in 
terms of the Hall coefficient £, the inverse plasma beta and the free parameter of the model go 
(i.e. the parameters which are expressed through the input characteristics values of the equilibrium 
disk, see Figs. 4, 5). 

The values of £ corresponding to the instability (see relations (29)) lie either between the 
branches £ _ and £ + for q > 1 (Fig. 4a, b for g > 0) or above the branch £~ for q < (Fig. 5a, b 
for g < 0). Then according to the data in Figs. 2 and 3, the disk is stable for = 0, but may 
become unstable for any finite plasma beta ft^ 1 > 0. Since ^ are inversely proportional to the 
derivative of the inertial moment density d r N m (see the second relation (28b) resolved with respect 
to £), the vertical asymptotes I in Figs. 4, 5 arise if d r N m =0 at the reference point r = 1. As seen 
in Figs. 4 and 5, both positive- and negative-^ modes are subdivided into two sub-modes 1 and 2 
separated by a vertical asymptote I, where d r N m =0, while the vertical asymptote II, where q = 0, 
bounds the region of instability (q < according to relations (29)) of the sub-mode 2. 

Figure 6 and estimations (24) demonstrate for typical values of the disk parameters that the 
Hall instability is excited with inverse growth rates of the order of the rotation period. Note that 
although the angular velocity of the Keplerian rotation drops out from the dimensionless stability 
problem for the toroidal magnetic configuration, the Keplerian rotation velocity determines the 
value of the Hall parameter (see Eqs. (6)). 

Let us estimate the minimal value of the dimensional toroidal magnetic field that corresponds 
to the onset of the Hall instability. According to the definition of £ in Eqs. (6), (20), the minimal 
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magnetic field that gives rise to the Hall instability, is given by (Fig. 7): 

4vren,(GM) 1 / 2 . 1/2 
-B* = ea^rj , (36) 

where £ may be estimated as the characteristic threshold value £~ that corresponds to the onset 
of the Hall instability in Figs. (4) and (5). Since in Eq. (36) is proportional to the ionization 
degree a and the aspect ratio of the disk e, the onset of the Hall instability may occur in ionized 
weakly and magnetized thin disks. 



5. Summary and conclusions 

In the present study, a linear analysis of the stratification driven Hall instability for thin 
equilibrium Keplerian disks with an embedded purely toroidal magnetic field and hydrostatically 
equilibrated density has been performed. The stability analysis has been carried out in the short- 
wave local approximation for the radial and vertical coordinates frozen at the reference point 
in the midplane, under the assumption of a small slope of the vertical edge of the disk. The 
leading order asymptotic expansions in small aspect ratio of the disk is employed in order to 
construct the Hall MHD equilibrium of the thin Keplerian disks. The Hall electric potential is 
introduced which determines the equilibrium toroidal magnetic field up to a free function of the 
inertial moment density. The latter has been modelled as a sum of the current free magnetic field 
(created by a current localized along the disk axis) and the simplest quadratic approximation for 
the magnetic field created by the electric current distributed over the Keplerial portion of the disk. 
The corresponding equilibrium profiles contain a single free parameter of the model that may be 
expressed through the total current - the natural parameter of the disk. 

It is shown that: (i) the radially-stratified equilibrium disk is stable to radial, and unstable 
under vertical perturbations. The instability is demonstrated by the marginal stability surface in 
the space of the local Hall and inverse plasma beta parameters, as well as the free parameter of 
the model (Fig. 1); (ii) the sign of the free parameter determines qualitatively different behavior 
of the both equilibrium and perturbations (Figs. 2-5). In particular, for negative values of the 
free parameter, the equilibrium disk exhibits a maximum in the number density that corresponds 
to rings of denser material, as well as a maximum in the magnetic field (Fig. 3); (hi) the disk is 
stable for the infinite value but may become unstable for any finite value of plasma beta (Figs. 2 
and 3); (iv) the current free configuration is shown to be stable, and the existence of the distributed 
electric current is found to be necessary in order to give rise to the Hall instability; (v) the inverse 
growth rate is of the order of the rotation period taken as the characteristic time scale (Fig. 6); 
(vi) the density inhomogeneity length is of the order of the characteristic disk thickness; (vii) the 
onset of the Hall instability is possible in thin weakly magnetized and ionized disks (Fig. 7). 

This work has been supported by the Israel Science Foundation under Contract No. 265/00. 
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Fig. 1. — The upper/bottom branches ( dotted/solid lines X + /X~) of the marginal stability curves 
in the plane of the local Hall and inverse plasma beta parameters b~ l — X for the reference value 
of radial coordinate r = 1 in the midplane z = (curves I, II, III for three different values g = go). 
Instability region lies either (a) between the upper and bottom branches {g > 1 and q > 1) or 
(b) above the bottom branch (g < and q < 0), dashed straight lines are the vertical asymptotes 
q = 0. 
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Fig. 2. — Equilibrium midplane profiles of the toroidal magnetic field (a), number density (b), 
distributed electric current (c) and inertial moment density (d) vs r for g = go = 3.5 and three 
values of the plasma beta (3 (curves I, II, III); £~ < £ = < C + belongs to the region of instability. 
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Fig. 3. — Equilibrium midplane profiles of the toroidal magnetic field (a), number density (b), 
distributed electric current (c) and inertial moment density (d) vs r for g = go = —0.8 and three 
values of the plasma beta (3 (curves I, II, III); £~ < £ = £o belongs to the region of instability. 
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Fig. 4. — Marginal stability curves in the plane of the Hall criterion and inverse plasma beta 
[3~ l — £ for (a) g = 3.5 and (b) g = 6 (g = go > 1, £ = £o, reference radius r = 1 in the midplane 
z = 0). Dotted/solid lines depict the upper/bottom branches £ = ^ of the marginal stability 
curves. Dashed lines depict the asymptote I (d r N m = 0) that separates the sub-modes 1 and 2. 
The upper branch in Fig. 4b and both branches in Fig. 4a for the sub-mode 1 are not presented 
since they correspond to the Hall parameter £ > 1.5. 
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Fig. 5. — Marginal stability curves in the plane of the Hall criterion and inverse plasma beta 
/3 _1 — £ for (a) g = —0.8 and (b) g = —1.2 ( g = go < 0, £ = £o, reference radius r = 1 in 
the midplane z = 0). Solid lines depict the marginal stability curves. Dashed lines depict the 
asymptote I (d r N m = 0) that separate the sub-modes 1 and 2, and the asymptote II (q = 0) which 
restricts the region of instability for the sub- modes 2. 
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Fig. 7. — Dimensional characteristic toroidal magnetic field in milli-Gauss vs reference radius in 
astronomical units for three values of the ionization degree a; £ = £o = 0.2; e = 0.02 ; = n n = 
10 12 cm- 3 ; GM = 1.3 • 10- 12 cm 3 /s 2 . 



